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Abstract. We have studied the muon neutrino and antineutrino quasi-elastic (QEL) scattering reac- 
tions {v^n — » n^p and i>^p — > j-i^n) using a set of experimental data collected by the NOMAD col- 
laboration. We have performed measurements of the cross-section of these processes on a nuclear tar- 
get (mainly Carbon) normalizing it to the total i/^ (i/^) charged current cross-section. The results for 
the flux averaged QEL cross-sections in the (anti)neutrino energy interval 3 — 100 GeV are {aqei)^,,. = 
(0.92±0.02(sfat)±0.06(sysf)) x 10"^* cm^ and = {0.81±0.05{stat) ±0.08{syst)) x 10"^* cm^for 

neutrino and antineutrino, respectively. The axial mass parameter Ma was extracted from the measured 
quasi-elastic neutrino cross-section. The corresponding result is Ma = 1-05 ± 0.02{stat) ± 0.06{syst) GeV. 
It is consistent with the axial mass values recalculated from the antineutrino cross-section and extracted 
from the pure shape analysis of the high purity sample of quasi-elastic 2-track events, but has 
smaller systematic error and should be quoted as the main result of this work. Our measured Ma is found 
to be in good agreement with the world average value obtained in previous deuterium filled bubble cham- 
ber experiments. The NOMAD measurement of Ma is lower than those recently published by K2K and 
MiniBooNE collaborations. However, within the large errors quoted by these experiments on Ma, these 
results are compatible with the more precise NOMAD value. 

Key words, neutrino interactions, neutrino and antineutrino quasi-elastic scattering 



PACS. 13.15.+g Neutrino interactions - 25.30.Pt Neutrino-induced reactions 
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1 Introduction 

A precise knowledge of the cross-section of (anti)neutrino- 
nucleus quasi-elastic scattering process (QEL) is impor- 
tant for the planning and analysis of any experiment which 
detects astrophysical, atmospheric or accelerator neutri- 
nos. The available measurements from early experiments 
at ANL miSpp], BNL [51I61I71I8] . FNAL [9l[H], CERN [111 
[HlfTSlfMlflKllTHlfTTlfTS] and IHEP[T9l[2(H[2T|[22] have con- 
siderable errors due to low statistics and a lack of knowl- 
edge of the precise incoming neutrino flux. Unfortunately, 
even within these large errors, the results are often con- 
flicting. 

This subject remains very topical. Recently several at- 
tempts have been made to investigate the QEL process 
in the data collected by modern accelerator neutrino ex- 
periments (such as NuTeV [23], K2K [241125] and Mini- 
BooNE [2S]). Unfortunately they have not clarified the sit- 
uation again due to large errors assigned to their measure- 
ments. Dedicated experiments, such as e.g. SciBooNE [27] 
and MINERi^A [5S], are now being performed. 

In the present analysis, we study both i^^ and QEL 
scattering in the data collected by the NOMAD collab- 
oration. The NOMAD detector was exposed to a wide- 
band neutrino beam produced by the 450 GcV proton 
synchrotron (SPS, CERN). A detailed description of the 
experimental set-up can be found in [55]. The character- 
istics of the incoming neutrino flux are given in [30| . 

The large amount of collected data and the good qual- 
ity of event reconstruction in the NOMAD detector pro- 
vide a unique possibility to measure the QEL cross-section 
with a combination of small statistical and systematic er- 
rors. The data sample used in this analysis consists of 
about 751000(23000) t-^ (Df,) charged-current (CC) inter- 
actions in a reduced detector fiducial volume. The average 
energy of the incoming {v^j) is 25.9 (17.6) GeV. 



The merit of the current analysis is the possibility of 
keeping the systematic error relatively small. It takes ad- 
vantage from three main factors: 1) the NOMAD detector 
is capable of selecting a sample of QEL events with a high 
purity and a good efficiency; 2) a simultaneous measure- 
ment of both two track and single track z^^ QEL events 
allows to constrain the systematics associated with nuclear 
reinteractions; 3) a wide energy range of the NOMAD neu- 
trino beam allows to perform a precise normalization to 
the well-known total (DIS) v^^ CC cross-section. 

The paper is organized as follows. In Section [2] we give 
a brief review of the published experimental data on QEL 
(anti)neutrino scattering. The NOMAD detector and the 
incoming neutrino flux are briefly discussed in Section [31 
In Section [4] we outline the MC modeling of signal and 
background events, emphasizing also the importance of 
nuclear effects. Section [5] is devoted to the selection of the 
QEL events; we describe the QEL identification proce- 
dure and compare the MC predictions with experimental 
data. The methods used to measure the QEL cross-section 
and the phenomenological axial mass parameter Ma are 
the subjects of Section[51 The systematic uncertainties are 
summarized in Section [T] The results are presented in Sec- 
tion [8j Finally, a summary and discussion of the obtained 
results are given in Section [S] 

More detailed information can be found in [31] . 



Correspondence to: Vladimir. Lyubiishkin@cern.ch 
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2 Review of existing experimental data on 
Quasi-Elastic (anti)neutrino scattering 

Let us start with a brief review of existing experimental 
data on (anti)ncutrino nucleon QEL scattering. 

A compilation of available data on the cross-section 
measurement of the z^^ and Dfj, quasi-elastic scattering 
off deuterons and other nuclei or composite targets (like 
freon, propane, liquid scintillator) as a function of the in- 
coming neutrino energy has been made (see Figures [T5 l II 61 
and fT7|) . This study allowed to conclude that the QEL 
cross-section measured in different experiments can vary 
by 20-40%. 

The existing data on (anti)neutrino QEL scattering 
come mostly from bubble chamber (BC) experiments. In 
general, these data suffer from small statistics. Moreover, 
results of several old experiments [T2[|13[|T4] have large 
systematic uncertainties due to the poor knowledge of the 
incoming neutrino flux and of background contamination 
in the selected events. 

The total QEL cross-section was recently measured in 
data collected by the NuTeV collaboration [23] . The num- 
ber of QEL events identified in their analysis are compara- 
ble with the total world data obtained in previous experi- 
ments. However, the results reported for the antineutrino 
case fall well outside the most probable range of values 
known today and hence, seem to exhibit a systematic shift. 

Another intriguing subject in the study of the neu- 
trino quasi-elastic scattering is the axial structure of the 
nucleon. We will skip here the details of the phenomenol- 
ogy of the hadronic current involved in the matrix element 
of the process (see Section |4T1 and Ref. [32]). But let us 
only remind the reader that for the region of low and in- 
termediate 4- momentum transfer, Q^, we can use a dipole 
parametrization for the axial form factor with only one 
adjustable parameter, the so-called axial mass Ma- 

The Ma parameter describes the internal structure of 
the nucleon and should be the same both for neutrino 
and antineutrino experiments (if we assume the isotopic 
invariance of strong interaction). Therefore, it is conve- 
nient to compare experimental results in terms of the ax- 
ial mass. There is, however, no theoretical basis for this 
form of the axial form factor. The use of an inappropriate 
parametrization could therefore lead to values of Ma that 
differ when extracted under different kinematical condi- 
tions. 

There are two possible ways generally used to extract 
the Ma parameter from experimental data: 

1. from the total QEL (anti)ncutrino nucleon cross-section 
(the axial form factor is responsible for about 50-60% 
of the total QEL cross-section); 

2. from the fit of the distribution of the identified 
neutrino QEL events. 

In principle, these two procedures should give self- 
consistent results. However, the old bubble chamber ex- 
periments at ANL and CERN reported in general larger 
values of Ma based on the fit than those obtained from 
the total cross-section measurements. 



Results of the Ma measurements based on the fit 
have been recently published by the K2K [241125] and Mini- 
BooNE [26] collaborations. They are about 15% higher 
than the average of previous deuterium filled bubble cham- 
ber experiments. This disagreement is, however, just at 
about one sigma level because of the large systematic er- 
rors associated with the K2K and MiniBooNE measure- 
ments. 

Let us note that the extraction of from the 
distribution fit is a more delicate issue than the QEL total 
cross section measurement. 

In general, there are at least three aspects which can 
affect noticeably the Ma measurements: 

1. The nuclear effects can distort the expected distribu- 
tions of the measured kinematic variables (like the en- 
ergy of the outgoing nucleon). The neutrino- nucleus 
interactions should be described by a theoretical model 
suitable for the considered neutrino energy region. This 
is important both for MC modeling in present-day neu- 
trino experiments and for a proper interpretation of 
the results obtained earlier (with few exceptions for 
the deuterium filled bubble chambers). 

2. The correct determination of the background contam- 
ination from both deep inelastic scattering and single 
pion production in the selected events is important for 
experiments operating with intermediate and high en- 
ergy neutrino beams. 

3. The QEL reconstruction efficiency as a function of 
for two-track events is not expected to be a flat func- 
tion. It should drop both at small due to the loss 
of low energy protons and at large due to the loss 
of low energy muons. Effects which influence the ef- 
ficiency of the low momentum particle reconstruction 
should be carefully taken into account in the MC mod- 
eling of the detector response. 

Table [1] and [2] display the measured values of Ma from 
neutrino and antineutrino experiments (this compilation 
is also presented in graphical form in Fig. [T5)l . Whenever 
possible we provide also the Ma measured from the total 
cross-section. 

From the results described above one can conclude that 
the presently available experimental data on the neutrino 
QEL cross-section allow for a very wide spread of the 
axial mass values, roughly from 0.7 to 1.3 GeV. There- 
fore the reliability of a theoretical fit to these data is 
questionable and the uncertainty attributed to such a fit 
should go beyond the averaged experimental statistical 
accuracy. Nevertheless, the formal averaging of Ma val- 
ues from several early experiments was done by the au- 
thors of [33]: Ma = 1.026 ± 0.021 GeV. This result is 
also known as the axial mass world average value. Ac- 
cording to [M1I551I55] an updated world average value from 
I'^-Deuterium and pion electroproduction experiments is 
Ma = 1.014 ±0.014 GeV. 
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Table 1. A summary of existing experimental data: the axial mass AIa as measured in previous neutrino experiments. Numbers 
of observed events have been taken from the original papers; usually they are not corrected for efficiency and purity (the so-called 
QEL candidates). The axial mass value for the NuTeV experiment [23] was estimated from the published neutrino quasi-elastic 
cross section ((7^' = (0.94 ± Om{stat) ± 0.07{syst)) x 10"^** cm^); the systematic error for IHEP SKAT 90 [22] is 0.14 GeV. 



Experiment 



Target 



Events 



Method 



Ma, GeV Ref. 



ANL 69 


Steel 




da/dQ^ 
a 


1.05 ±0.20 
0.97 ±0.16 


m 


ANL 73 


Deuterium 


166 


da/dQ^ 
a ® da/dQ'^ 


0.94 ±0.18 










0.95 ±0.12 










a 


0.75±«;i? 




ANL 77 


Deuterium 


~ 600 


da/dQ^ 
a (g) da/dQ^ 


1.01 ±0.09 


m 








0.95 ±0.09 










a 


0.74 ±0.12 




ANL 82 


Deuterium 


1737 


dajdQ^ 
a ® dajdQ^ 


1.05 ±0.05 


m 








1.03 ±0.05 




BNL 81 


Deuterium 


1138 


dajdQ^ 
da/dQ^ 


1.07 ±0.06 


m 


BNL 90 


Deuterium 


2538 


1 070+°-°'*° 


m 



FermiLab 83 
NuTeV 04 
MimBooNE 07 



Deuterium 
Steel 

Mineral oil 



362 
21614 
193709 



da/dQ^ 
a 

da/dQ^ 



1.05j 



f0.12 
-0.16 

1.11 ±0.08 
1.23 ±0.20 



CERN HLBC 64 
CERN HLBC 67 
CERN SC 68 
CERN HLBC 69 



Freon 
Freon 
Steel 
Propane 



CERN GCM 77 Freon 

CERN GCM 79 Propane/Freon 

CERN BEBC 90 Deuterium 



236 
90 
236 
130 

687 

556 
552 



da/dQ^ 
a ® dajdQ"^ 

da/dQ^ 
a ® da/dQ^ 



i.oo: 



da/dQ^ 
a 

da/dQ^ 
a 

da/dQ^ 



fO.35 
-"-0.20 

''J-0.20 

U.DO_o.40 

0.70 ± 0.20 
0.88 ±0.19 
0.96 ±0.16 
0.87 ±0.18 
0.99 ±0.12 
0.94 ± 0.07 
1.08 ±0.08 



K2K 06, SeiFi Water 
K2K 08, SciBar Carbon 



12000 



da/dQ'^ 
da/dQ^ 



1.20 ±0.12 
1.144 ±0.077 



m 

[23] 

m 



m 

M 
M 

m 
m 

m 



IHEP 82 


Aluminium 


898 da/dQ'^ 
1753 da^+i,/dQ^ 
464 a®da/dQ' 
a 


1.00 ±0.07 


m 


IHEP 85 


Aluminium 


1.00 ±0.04 


[20] 


IHEP SCAT 88 


Freon 


0.96 ±0.15 
1.08 ±0.07 


m 


IHEP SCAT 90 


Freon 


da/dQ^ 
a (g) dajdQ^ 


1.05 ±0.07 


[22] 






1.06 ±0.05 





m 

[251 



Table 2. The same as in Table [T] but for antineutrino experiments. The axial mass value for the NuTeV experiment [23] was 
estimated from the published antineutrino quasi-elastic cross section (erf' = (1.12 ± 0.04(stat) ± 0.10(si/sf)) x 10"^* cm ); the 
systematic error for IHEP SKAT 90 [22] is 0.20 GeV. 



Experiment 


Target 


Events 


Determined from 


Ma, GeV 


Ref. 


BNL 80 
BNL 88 


Hydrogen 
Liquid scint. 


2919 


da/dQ^ 
da/dQ^ 


9+°* 
1.09 ±0.04 


m 
m 


FermiLab 84 
NuTeV 04 


Neon 
Steel 


405 
15054 


da/dQ^ 
a 


0.99 ±0.11 
1.29 ±0.11 


m 

[23] 


CERN GCM 77 
CERN GCM 79 


Freon 

Propane/Freon 


476 
766 


a 

da/dQ'^ 
a 

da/dQ'^ 


0.69 ±0.44 
0.94 ±0.17 

0.91 ±0.04 


m 
m 


IHEP 8.5 
IHEP SKAT 88 

IHEP SKAT 90 


Aluminium 
Freon 

Freon 


854 
52 


dov+v 1 dQ^ 
a g) dajdQ^ 
a 

dajdQ^ 
a (g) do/dQ"^ 


1.00 ± 0.04 
0.72 ±0.23 
0.62 ±0.16 
0.79 ±0.11 
0.71 ±0.10 


ED] 

m 

[22] 
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3 The NOMAD detector 

The NOMAD detector [29] consisted of an active target 
of 44 drift chambers with a total fiducial mass of 2.7 tons, 
located in a 0.4 Tesla dipole magnetic field as shown in 
Fig. [TJ The X xY x Z total volume of the drift chambers 
is about 300 x 300 x 400 cm^. 

Drift chambers [37j . made of low Z material served 
the dual role of a nearly isoscalar targelQ for neutrino in- 
teractions and of tracking medium. The average density 
of the drift chamber volume was 0.1 g/cm^. These cham- 
bers provided an overall efficiency for charged track re- 
construction of better than 95% and a momentum resolu- 
tion which can be approximated by the following formula 
^ « ^ © o^p)^ where the momentum p is in GeV/c 

and the track length L in m. Reconstructed tracks were 
used to determine the event topology (the assignment of 
tracks to vertices), to reconstruct the vertex position and 
the track parameters at each vertex and, finally, to iden- 
tify the vertex type (primary, secondary, etc.). A transi- 
tion radiation detector (TRD) [381139] placed at the end 
of the active target was used for particle identification. 
Two scintillation counter trigger planes [10] were used to 
select neutrino interactions in the NOMAD active target. 
A lead-glass electromagnetic calorimeter [4T1I42] located 
downstream of the tracking region provided an energy res- 
olution of 3.2%/ y/ E[GeV] © 1% for electromagnetic show- 
ers and was crucial to measure the total energy flow in 
neutrino interactions. In addition, an iron absorber and 
a set of muon chambers located after the electromagnetic 
calorimeter was used for muon identification, providing 
a muon detection efficiency of 97% for momenta greater 
than 5 GeV/c. 

The NOMAD neutrino beam consisted mainly of t'p's 
with an about 7% admixture of and less than 1% of 
Ve and De- More details on the beam composition can be 
found in [30] . 

The main goal of the NOMAD experiment was the 
search for neutrino oscillations in a wide band neutrino 
beam from the CERN SPS [43][44]. A very good quality 
of event reconstruction similar to that of bubble chamber 
experiments and a large data sample collected during four 
years of data taking (1995-1998) allow for detailed studies 
of neutrino interactions. 



3.1 Reconstruction of QEL events in the NOMAD 
detector 

A detailed information about the construction and perfor- 
mance of the NOMAD drift chambers as well as about the 
developed reconstruction algorithms is presented in [37]. 
Let us briefly describe some features relevant to the cur- 
rent QEL analysis. The muon track is in general easily 
reconstructed. However, when we study protons emitted 
in the Vfj, QEL two-track candidates we deal with protons 



Neutrino 



Dipole Magnet 
Front ® B = 0.4 T TRD 

Calorimeter | Modules Presliower 



nnnnn[|ir|n 




Muon 
Ciiambers 



Drift Chambers 



Electromagnetic Hadronic 
Calorimeter Calorimeter 



Fig. 1. A side-view of the NOMAD detector. 



with momentum well below 1 GeV/c and with emission 
angle above 60 degrees. For positive particles in the up- 
ward hemisphere of the NOMAD detector such conditions 
mean that these particles are almost immediately making 
a U-turn due to the magnetic field. There were no spe- 
cial efforts invested into tuning the NOMAD reconstruc- 
tion program to reconstruct this particular configuration 
(which is rather difficult due to the fact that these protons 
are in the l//3^ region of ionization losses, traversing much 
larger amount of material, crossing drift cells at very large 
angles where the spacial resolution of the drift chambers is 
considerably worse and where a large amount of multiple 
hits is produced, etc.). Some of these effects are difficult 
to parametrize and to simulate at the level of the detec- 
tor response in the MC simulation program. Thus, the 
reconstruction efficiencies for this particular configuration 
of outgoing protons could be different for the simulated 
events and real data. 

Let us stress, however, that for protons emitted down- 
wards we observed a good agreement between data and 
MC. 

In the current analysis it was important to disentangle 
the reconstruction efficiency effects discussed above from 
the effects induced by intranuclear cascade (which could 
change the proton kinematics and thus introduce drastic 
changes in the final results due to the efficiency mismatch 
between simulated and real data) . In order to get rid of an 
interplay between these two effects it was crucial to choose 
the region in the detector with a stable reconstruction effi- 
ciency. This could be achieved by selecting QEL events 
where protons are emitted in the lower hemisphere of the 
NOMAD detector. This approach allowed to find the best 
set of parameters for description of the intranuclear cas- 
cade. 

The most upsteam drift chamber was used as an addi- 
tional veto to remove through-going muons from neutrino 
interactions upstream of the NOMAD active target. This 
is crucial for the study of single track events. 



^ the NOMAD active target is nearly isoscalar (n„ : Up = 
47.56% : 52.43%) and consists mainly of Carbon; a detailed de- 
scription of the drift chamber composition can be found in [37] 
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4 Monte Carlo simulation of neutrino 
interactions 

Inclusive (anti) neutrino charged current (CC) and neutral 
current (NC) scattering can be considered as a mixture of 
several processes described by significantly different mod- 
els. In our case, these are quasi-elastic scattering (QEL), 
single-pion production (RES) and deep inelastic scatter- 
ing (DIS). There is also a contribution from a coherent 
interaction of neutrino with a nucleus as a whole (COH). 
Below we will describe in details the simulation scheme 
used for each of these processes and discuss the influence 
of the nuclear effects. 

An adequate MC description of neutrino interactions is 
important to calculate the efficiency of the QEL selection. 
Moreover, it allows us to predict the level of background, 
which cannot be suppressed completely by the QEL iden- 
tification scheme proposed in Section [5l 



4.1 Quasi-elastic neutrino scattering 

The standard representation of the weak hadronic current 
involved in the matrix elements of the processes i/^n — > 
H^p and V^p IJ^^n, is expressed in terms of 6 form- 
factors, which in general are assumed to be complex [451 . 
They formally describe the hadronic structure and can- 
not be calculated analytically within the framework of the 
electro- weak interaction theory. 

We neglect the second-class current contributions as- 
sociated with the scalar and pseudo-tensor form-factors. 
This is equivalent to the requirement of time reversal in- 
variance of the matrix element (hence all form-factors 
should be real functions of Q^) and charge symmetry of 
the hadronic current (rotation about the second axis in 
the isotopic space). 

The vector form-factors Fy and Fm are related through 
the isospin symmetry hypothesis to the electromagnetic 
ones, which we will consider to be well known. Instead of 
the simple dipole parametrization, extensively used in pre- 
vious experiments, we have chosen the Gari-Kriiempelmann 
(GK) model 05] extended and fine-tuned by Lomon [Tf] . 
Specifically we explore the "GKex(05)" set of parame- 
ters [48] which fits the modern and consistent older data 
well and meets the requirements of dispersion relations 
and of QCD at low and high 4-momentum transfer [46] . 

For the axial and pscudoscalar form factors we use the 
conventional representations [45j : 



spread. Thus, in our case, it should be considered as one 
of the available parameters, which can be used to adjust 
the MC simulation with the measured value of the total 
QEL cross section and observed distributions of the kine- 
matic variables (other parameters, related to the modeling 
of the intranuclear cascade, will be described later). 

Note that the expression for the pscudoscalar form fac- 
tor Fp is nothing better than a plausible parametrization 
inspired by the PCAC hypothesis and the assumption that 
the pion pole dominates at ^ 'tIt S5]- However, its 
contribution enters into the cross sections multiplied by 
a factor (m^/mjv)^. Hence, the importance of the related 
uncertainty is much reduced. 



4.2 Single-pion production through intermediate 
baryon resonances 

In order to describe the single-pion neutrino production 
through baryon resonances we adopt an extended version 
of the Rein and Sehgal model (RS) [50ll51j . which seems 
to be one of the most widely trusted phenomenological 
approaches for calculating the RES cross sections. The 
generalization proposed in ^52.. 53] takes into account the 
final lepton mass and is based upon a covariant form of 
the charged leptonic current with definite lepton helicity. 
In our MC simulation we use the same set of 18 interfering 
nuclcon resonances with masses below 2 GeV as in [50] but 
with all relevant input parameters updated according to 
the current data [491154] . Significant factors (normalization 
coefficients etc.), estimated in Ref. [SO] numerically are 
recalculated by using the new data and a more accurate 
integration algorithm. 

The relativistic quark model of Feynman, Kislinger, 
and Ravndal [55j . adopted in the RS approach, unam- 
biguously determines the structure of the transition ampli- 
tudes involved into the calculation and the only unknown 
structures are the vector and axial- vector transition form 
factors G^'^(Q2). In [50] they are assumed to have the 
form 



and 



Fa{Q^)^Fa{Q)[1 



Fp (g^) = 



2m^ 



Ml 



(1) 



(2) 



where Fa{0) = gA = —1.2695 ± 0.0029 (measured in neu- 
tron /3-decay |49| ) : and m^y - pion and nuclcon masses. 

As discussed in Section [21 the currently available ex- 
perimental data on the axial mass AIa allow for a wide 



4m% 



l/2-n 



M2 



(3) 



where the integer n in the first ( "ad hoc" ) factor in Eq. ([3]) 
is the number of oscillation quanta of the intermediate 
resonance. 

The vector mass My is taken to be 0.84 GcV, that is 
the same as in the usual dipole parametrization of the nu- 
clcon electromagnetic form-factors. The axial mass (which 
was fixed at 0.95 GcV in the original RS paper) is set to 
the standard world averaged value AIa = 1-03 GeV. It is 
in good agreement with the results obtained in the recent 
analysis of the data from the BNL 7-foot deuterium filled 
bubble chamber [5^ {Ma = 1.08 ±0.07 GeV). Let us also 
note that the available experimental data for the single- 
pion neutrino production (as in the case of QEL scatter- 
ing) does not permit a very definite conclusion about the 
value of the total RES cross section (and the correspond- 
ing axial mass value). The present uncertainties will be 
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taken into account in the calculation of the systematic 
error of the current analysis. 

To compensate for the difference between the SUq pre- 
dicted value (—5/3) and the experimental value for the 
nuclcon axial- vector coupling gA, Rein and Sehgal intro- 
duced a renormalization factor Z — 0.75. In order to ad- 
just the renormalization to the current world averaged 
value gA = -1.2695 [49] we have adopted Z 0.762. 
The harmonic-oscillator constant fi, which accounts for 
the mass differences between states with different num- 
ber of excitation quanta is set to its original value Q = 
1.05 GcV^ 

Another essential ingredient of the RS approach is 
the non- resonant background (NRB). Its contribution is 
important in describing the existing data on the reac- 



tions I'fj^n 



fi riTT"*", ly^n — > /i pir'^ , v^p ^^pir and 



i^fiP — *■ fj, ' n7r°. In our Monte Carlo, the NRB is taken to 
come from the DIS part of the simulation. Therefore it 
has not been used in the RES part of our event generator. 
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4.3 Deep inelastic scattering 



Fig. 3. Benhar-Fantoni parametrization [70] for the momen- 
tum distribution of the target nucleons (solid line) , normalized 
to the Fermi distribution with zero temperature and Fermi 
momentum Pp = 221 MeV/c (simple RFG, dashed line). 



The MC simulation of the deep inelastic neutrino nuclcon 
scattering is based on the LEPTO 6.5.1 package [S^ with 
several modifications [59ll60] . For hadronization we use the 
LUND string fragmentation model, as incorporated into 
the JETSET 7.4 program [^[MIIMI- 

Upon implementing the Monte Carlo for v^{v^) CC 
scattering, kinematic boundaries between exclusive (RES) 
and inclusive (DIS) channels must be defined. To avoid 
double counting, the phase space of the RES and DIS 
contributions should be separated by the conditions W < 
W^f^ and W > W^(^ , where W is the invariant mass 
of the final hadronic system. 

The maximum possible value for W^f^ is the upper 
limit of the RS model (2 GeV), while inelastic scatter- 
ing can take place from the one-pion production thresh- 
old (note, however, that this value is too small in principle 
since the structure functions used in the calculation of the 
DIS cross section cannot be extrapolated down to this 
value) . 

Unfortunately, there is no clear physical recipe to de- 
termine exact numerical values for those cutoff parame- 
ters. The authors of GENIE MC code [Bi] adopt the value 
W^f^ ~ W^'^ ~ 1.7 GcV. A comprehensive analysis of 



4.4 Coherent pion production 

In the processes described above, neutrinos interact with 
individual target nucleons. However, pions can be pro- 
duced in a coherent interaction of the neutrino with the 
whole nucleus, i.e. in the case of CC i^,, scattering M 
^~TT'^N , where N is the target nucleus. 

The details of the MC simulation can be found in [SB] , 
which is devoted to the investigation of this process in the 
NOMAD experiment. The flux averaged cross-section has 
been calculated following [67l[68] and has been estimated 
at 0.733 X lO^'^^cm^ per nucleus. For a recent experimen- 
tal result at low incoming neutrino energy see [69| . Taking 
into account that the average mass number of the NO- 
MAD target is 12.9, and using the number of recorded 
DIS events (see section [B.l.ip one finds that the expected 
number of coherent pion production events is ~ 2700. Nev- 
ertheless, the probability for events of this type to be iden- 
tified as QEL is ~ 2% because of the small pion emission 
angle, so that the expected contamination of the selected 
QEL sample is lower than 0.4%. 



available experimental data made in [57|l65] suggests to 
decrease this cut to ~ 1.5 GeV. 

In the present analysis we set W^^^ = 2 GeV and 
Wl^Jt^ = 1.4 GeV. This choice allows for the non-resonant 
contribution to single pion production to be accounted 
for by the DIS part of the Monte Carlo, providing e.g. 
N{fj.-pTr° in mS)/Ndis-<Jdis » (j{fi-pTr° in NRB from RES) 
see previous subsection. Moreover, it is not at variance 
with experimental data as far as the total (anti)neutrino 
cross-section is concerned (see Fig. [2]). 



4.5 Nuclear effects 

For typical NOMAD neutrino energies, we can assume 
that the incident neutrino interacts with one nucleon only 
inside the target nucleus, while the remaining nucleons 
,are spectators (Impulse Approximation). In this case, one 
can describe the neutrino nucleus scattering by folding 
the usual expressions for the free neutrino nucleon cross 
sections with a Fermi gas distribution. 
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Fig. 2. Slopes of the total i/^ and CC scattering cross-sections off an isoscalar nucleon (the compilation of experimental 
data is taken from [57])- The curves and bands show the QEL, RES, and DIS contributions and their sums calculated with the 
parameters described in the legend of the top panel. The averaged values over all energies (0.677 ± 0.014) x 10^^® cm^/GeV (for 
iy^,N) and (0.334 ± 0.008) x 10"^** cm^/GeV (for V^N) obtained by the Particle Data Group [49] are also shown for comparison 
(straight lines). 



9 



In the rclativistic Fermi gas model, the nucleus is con- 
sidered as an infinite system of non-interacting nucleons. 
The phenomena related to the nuclear surface and to the 
interaction between nucleons can be taken into account by 
using a more realistic effective momentum distribution for 
the target nucleons. In the NOMAD event generator we 
used the Benhar-Fantoni parametrization [70], see Fig.[3l 

The QEL simulation is based on the Smith-Moniz ap- 
proach [71| . The momentum of the recoil nucleus and the 
nucleon binding energy are included in the conservation 
laws which determine the event kinematics. The only final 
state interaction (FSI) effect which is taken into account 
at this stage is the Pauli exclusion principle. The explicit 
form of the QEL differential cross section used in the MC 
code can be foimd in [32] . 

MC implementation of the Fermi gas model in the case 
of single pion production is more straightforward. First, we 
generate the momentum of the target nucleon and make 
a Lorentz boost to its rest frame where the RES event 
can be simulated according to the extended RS model de- 
scribed in subsection 14.21 The effect of Pauli blocking on 
the outgoing nucleon is taken into account as it is in the 
QEL MC. 

In the case of the DIS neutrino scattering there are 
several specific nuclear effects (such as nuclear shadow- 
ing, pion excess and off-shell corrections to bound nucleon 
structure functions) . They are described in the theoretical 
framework proposed in [72j . 

Simulating the re-interactions between particles pro- 
duced at the primary neutrino collision off the target nu- 
cleon with the residual nucleus is an important ingredient 
of the MC event generator. To include this effect, com- 
monly called final state interactions, we use the DPMJET 
package [73] . 

The intranuclear re-interaction of the particles gener- 
ated by the QEL, RES or DIS event generators can be de- 
scribed and simulated by the Formation Zone Intranuclear 
Cascade model [74 1 175] implemented in DPMJET. Secon- 
daries from the first collision are followed along straight 
trajectories and may induce in turn intranuclear cascade 
processes if they reach the end of their "formation zone" 
inside the target; otherwise they leave the nucleus without 
interacting. 

There are two important parameters in DPMJET. The 
first one, called the formation time tq, controls the de- 
velopment of the intranuclear cascade. With increasing 
To, the number of cascade generations and the number of 
low-energy particles will be reduced. Its default value is 
To = 2.0. After some tuning described below we adopted 
the value tq = 1.0 in our simulation of QEL, RES and DIS 
events. 

Inside DPMJET, the momenta of the spectator nucle- 
ons are sampled from the zero temperature Fermi-distribu- 
tion. However, the nuclear surface effects and the interac- 
tion between nucleons result in a reduction of the Fermi 
momentum, see Fig. |3J It can be accounted for by intro- 
ducing a correction factor a^^^^ (default value 0.6). More- 
over, a£!„jj provides the possibility of some modification 
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Fig. 4. Flux averaged cross-section of QEL (anti) neutrino 
scattering for NOMAD v^iip^i) beam as a function of the axial 
mass Ma. 



of the momentum distribution for the emitted low-energy 
nucleons. 

At the end of the intranuclear cascade, the residual nu- 
cleus is supposed to go through some de-excitation mech- 
anisms. It can be disrupted into two or more fragments, 
emit photons, nucleons or light particles (like d, a, '^H, 
■^He). We can easily neglect this contribution, since the 
typical energy of those particles is below the observation 
threshold of the NOMAD detector. 

In our analysis, special attention will be devoted to the 
dependence of the obtained results on the intranuclear cas- 
cade parameters. As a cross-check, we compare our MC 
simulation for the QEL process with the predictions of the 
NUANCE event generator [76], which is currently used in 
a large number of neutrino experiments and which con- 
tains a different approach to the modeling of FSI effects. 



4.6 Expected signal/background ratio in the 
^^(^A') sample 

In this subsection we estimate the number of signal quasi- 
elastic events in the initial i'fj,{v^) CC sample. 

The contribution of each process to the total set of 
events is proportional to its flux averaged cross-section: 

(a) = J G{E,)<P{E,)dE, I J ^{E,)dE, (4) 



where 



{E 



is the theoretical prediction for the cross-section of the 
process at stake, (1>{E^) denotes the NOMAD (anti)ncutrino 
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Table 3. Flux averaged cross-sections of the QEL, RES, DIS 
CC and NC processes per one nucleon of the NOMAD target. 
Neutrino beam spectrum corresponds to the |X, y| ^ 100 cm 
fiducial area. The unit used for the cross-section is 10~^* cm^. 



Process type 






QEL 


0.428 


0.393 


RES 


0.576 


0.432 


DIS CC 


16.643 


4.876 


DIS NC 


5.335 





energy spectrun:|3; nn{np) is the relative fraction of neu- 
trons (protons) in the NOMAD target (see Section [3]). 

The QEL cross-section was calculated in the frame- 
work of the Smith and Moniz model |71j for Carbon with 
binding energy — 25.6 McV and Fermi momentum 
Pp = 221 MeV/c. As noted above, the final result de- 
pends strongly on the axial mass Ma (see Fig.|4]). 

To estimate the RES contribution, wc fold the ex- 
tended RS model [53] for a free nucleon with the Pauli 
factor from [77]. The computation of (Jdis{Eu) has been 
done with the GRV98-LO PDF model as indicated in [57] , 
The cutoff parameters W^^^ and W^l^l^ are the same as 
for the MC simulation. 

Table [3] contains our results for the reduced fiducial 
volume of the NOMAD detector: \X, Y\ ^ 100 cm; the 
average (i/^) energy was 25.9 (17.6) GeV. 

Combining all these, the expected fraction of quasi- 
elastic events in the initial Vp,{y^) CC sample before any 
special selection is about 2.4%(6.9%) or - 2G3GG(- 1360) 
events. 



5 Events selection 

In this section we describe particular features of recon- 
struction and identification of and QEL events. 

5.1 u^n jjL^p selection 

For a v^n — > pTp event one can expect two tracks originat- 
ing from the reconstructed primary verteslfj: one of them 
should be identified as a muon, while the second track is 
assumed to be a proton. Later we shall refer to events with 
such a topology as 2-track (two track) eventfl 

Sometimes the proton track cannot be reconstructed, 
e.g. if its momentum is below the detector registration 

^ the procedure used for the calculation of the flux and com- 
position of the CERN SPS neutrino beam is described in [30] 

^ all charged tracks originating within a 5 cm box around the 
reconstructed primary vertex are forced to be included into it; 
we have also tried to vary this parameter by enlarging the size 
of the box to 10 cm and found that the final results are rather 
stable (within 0.3% for the measured QEL cross-section) 

* in this analysis we do not take into account clusters in 
the electromagnetic calorimeter, which can be associated with 
neutral particles, originating from the primary vertex 



threshold. In this case, we deal with only one muon track 
and we call such an event a 1-track (single track) event. 

The expected ratio between 1-track and 2-track events 
for the pure standard QEL MC sample is 54.3% : 45.7%. 

There are three possible reasons for the reconstruction 
of the proton track in a QEL event to fail: 

— the proton, which was born in the neutrino interaction 
with the target nucleon, has too low a momentum or 
too large an emission angle (this depends on the pa- 
rameters of the model used to describe the neutrino- 
nucleon interaction, in particular, on the value of the 
axial mass); 

— the proton from the primary neutrino interaction was 
involved in an intranuclear cascade and lost part of its 
energy (this is controlled by the DPMJET parameters, 
mainly by the formation time rg); 

— the detector magnetic field deviates positively charged 
particles upwards; therefore, if a slow proton is emit- 
ted at an azimuth iph ~ 7''/2, its trajectory is almost 
parallel to the drift chamber planes and its track re- 
construction efficiency (which depends on the number 
of hits associated with the track) is significantly lower 
than in the case of a proton emitted downwards at 
iph ^ 37r/2. 

In Fig. [5] (left) we illustrate these last two effects: the 
magnetic field is the cause of the asymmetry in the azimu- 
tal distribution of the reconstructed protons, while variyng 
the formation time parameter tq affects the expected num- 
ber of tracks uniformly. 

In Fig. [6] we display an example of distributions of the 
leading proton momentum ph and emission angle 9h before 
and after FSI for the QEL neutrino scattering. The proton 
reconstruction probabilities are also shown as functions of 
Ph and Oh'- one can observe a fast decrease at low proton 
momenta (below 300 MeV/c) and large emission angles 
(larger than 72°). So, FSI tends to increase the fraction 
of events in kinematic domains with low proton recon- 
struction efficiency and therefore to change the expected 
fraction of events with a given topology in the identified 
QEL sample. 

Using 2-track events only for the analysis may seem 
very attractive, since we could significantly reduce the 
background contamination with the help of additional kine- 
matic variables (details can be found below). However, 
the results thus obtained might still have large system- 
atic uncertainties coming from insufficient understanding 
of nuclear effects. 

The QEL events which are not reconstructed as 2- 
track events will populate mainly the 1-track sample. But 
(jqei extracted from this sample will suffer from the same 
source of uncertainty. However, the measurement of the 
QEL cross-section simultaneously from both samples is 
expected to have only little dependence on the uncertain- 
ties in the modeling of FSI effects and this is indeed what 
is found in the data (see Section [5]). 

Therefore, the strategy of our analysis (selection crite- 
ria) in the case of Vf^n fi^p can be outlined as follows: 
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Fig. 5. The reconstruction efficiency of proton track as a function of its azimuth tpu for QEL scattering; the curves are 
smoothed MC predictions obtained for different values of the formation time ro (left). The comparison of the muon azimuth (/^^ 
distributions in data and MC for 1-track and 2-track samples (right). 
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Fig. 6. Distribution of the leading proton momentum (left) and emission angle (right) before (dash-dotted line) and after 
(solid line) FSI simulation. Dashed lines show the proton reconstruction efficiency as function of the proton momentum and 
emission angle (for ir < iph < 27r). 



— Fiducial volume cut. The reconstructed primary vertex 
should be within the restricted fiducial volume (FV): 

\X, Y\ s$ 100 cm, 25 Z ^ 395 cm (5) 

— Identified muon. We require the presence of a recon- 
structed and identified negatively or positively charged 

^ we use a more stringent cut Z > 50 cm for the data col- 
lected during 97 and 98, when the first drift chamber module 
was substituted by the NOMAD STAR detector 



muon for the neutrino and antincutrino analyses re- 
spectively. In order to avoid possible problems with 
detector reconstruction inefficiencies, we require < 
iPf^ < TT, where (p^ is the muon azimuthal angle (so, 
the proton track should lie in the bottom hemisphere), 
see Fig. [5] (right). 

This choice is validated by our final errors being dom- 
inated by systematics as we will be shown below. 



12 



Table 4. Number of data Ndata and renormalized MC Nmc events in i/^ and QEL samples; expected selection efficiency, 
purity and background contaminations (BG) for different stages of the analysis. 







en. 


lo^rj-L purity 


rlH/o ovjr \m) 


JJlo JjLt 


Otn. oLj 


AT 
' data 


AT 
' mc 




1-track before Qh cut 


23.7 


29.0 


18.3 


52.3 


0.4 


16508 


16633.7 




1-track after 6^^ cut 


21.3 


41.7 


23.2 


34.5 


0.6 


10358 


10358.0 




2-track before L cut 


17.6 


47.2 


17.3 


35.2 


0.3 


7575 


7609.0 




2-track after L cut 


13.3 


73.9 


10.2 


15.8 


< 0.1 


3663 


3663.0 




combined before cuts 


41.3 


34.7 


18.0 


47.0 


0.3 


24083 


24242.7 




combined after cuts 


34.6 


50.0 


19.8 


29.7 


0.5 


14021 


14021.0 




1-track before cut 


81.8 


29.8 


22.8 


45.8 


1.6 


3585 


3555.8 




1-track after Qh cut 


64.4 


36.6 


28.5 


33.6 


1.3 


2237 


2237.0 



— Event topology and reconstructed kinematic variables. 
We assign the events to the 1-track and 2-track sub- 
samples and calculate E^, and Q^. 

• Single track sample (only one charged lepton is re- 
constructed and identified). To avoid contamina- 
tion from the through-going muons we extrapolate 
the muon track to the first drift chamber and re- 
quire the absence of veto chamber hits in the vicin- 
ity of the intersection point. The efficiency of this 
quality cut was controlled by visual scanning of 
the reconstructed 1-track events in the experimen- 
tal data and was found to be satisfactory. Another 
quality cut was used to suppress a possible contri- 
bution from inverse muon decay events: we require 
the muon transverse momentum to be greater than 
0.2 GeV/c (see Section [01 for more details). 

The kinematic variables are reconstructed under 
the assumption that the target nucleon is at rest. 
For the 1-track events, the muon momentum and 
direction are the sole measurements and we have 
to use the conservation laws (assuming QEL) to 
compute other kinematic quantities: 

^ ME^ - ml/2 

M - Efj^+ pi_, cos Ofj, 
g2 = 2MiE, - E^) 

Ph - ((S. -p^COS0,,)2 gjj^2 Q^Y'^ 

cos 6h = {Ey - p^ cos e^)/ph, (6) 

where p^, 9^ (p/j, Oh) are the momentum and emis- 
sion angle of the outgoing muon (nucleon), see Fig.[71 
We note that for the neutrino energies relevant 
for this analysis (above 3 GeV) there is no dif- 
ference between the calculations based on the ap- 
proximated formulae above and the precise one, 
which takes into account the binding energy (see 
e.g. Eq. (4) in [26]). With the help of the MC 
simulation we estimate the resolution of the recon- 
structed El, and as 3.6% and 7.8% respectively. 

• Two track sample (both the negative muon and the 
positively charged track are reconstructed). For a 
reliable reconstruction, we require that the number 
of hits associated with the positively charged track 




Fig. 7. Likelihood variables: missing transverse momentum 
P]"", proton emission angle 9h, angle a between the transverse 
components of the charged tracks. 

should be greater than 7 and its momentum ph > 
300 MeV/c. Otherwise such an event is downgraded 
to the 1-track sample. 

For 2-track events, we use both the muon and the 
proton reconstructed momenta to estimate E^, and 

El, = Pf^ cos 0f_i + Ph cos Oh 

Q2 ^ 2E,{E^ - p^ cos 0^)-ml 

The expected resolutions for Ei, and are 3.6% 
and 7.1%. 

The quality of the neutrino energy Ei, reconstruction 
for 1- and 2-track samples is illustrated in Fig.[8l It was 
checked that for the 2-track sample the derived cross- 
sections are consistent within errors for both methods 
of El, calculation. 
— Background suppression. The contamination from RES 
and DIS processes can be suppressed by using the dif- 
ference between kinematical distributions in the QEL 
and background events as well as by the identification 
of the reconstructed positively charged track as a pro- 
ton (for the 2-track sample only) . Therefore we apply: 
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Fig. 8. The quality of the neutrino energy Ei, reconstruction for 1- and 2-track samples. 



• Identification of the positively charged track. 
Momentum-range method |78] can be reliably ap- 
plied for low energy protons since their tracks are 
shorter compared to that of 7r+ (the main back- 
ground for proton identification) due to larger ion- 
ization losses. In our case, this method can be ap- 
plied to about 17% of the events El- 

• Kinematical criteria. 

In the case of the 2-track sample, we can use addi- 
tional kinematic variables to suppress background 
contamination. We build the likelihood ratio 



£ = ln 



P{1 \QEL) 
P{1 \BG) ' 



(7) 



using 3-dimensional correlations between the fol- 
lowing kinematic variables (see Fig. [7]): the missing 
transverse momentum p™"^ the proton emission 
angle Oh, and the angle a between the transverse 
components of the charged primary tracks. The fol- 
lowing pre-cuts were applied prior to the likelihood 
construction: P™'' < 0.8 GeV/c, 0.2 Ou/t^ ^ 0.5 
and a/iT 0.8. 

In Eq. the P(£ \QEL) and P(£ \BG) are the 
probabilities for signal and background events to 
have the values of the variables £. = (P™**, 6*/!, a). 



® We also undertook an attempt to identify positively 
charged particles using the TRD information. A special al- 
gorithm [79ll80j can be potentially used for discrimination be- 
tween two particle-ID hypotheses (p/tt in our case). However, 
a low momentum (~ 0.9 GeV) of the particle and a rather 
large emission angle (> 45°) result in that either the particle 
does not reach the TRD or the number of residual TRD hits 
is not large enough for the identification. Therefore, the TRD 
algorithm could be applied only to a limited fraction of events 
(~ 6%) and cannot play any significant role in our analysis. 



We have found that the DIS and RES probability 
functions are very similar; therefore we build the 
likelihood function taking only resonance events for 
the denominator of Eq. ([7|) . 

The comparison of P™''', a, 9h and C distributions 
in the data with the proper mixture of simulated 
QEL, RES and DIS events is displayed in Fig. [HI 
The good agreement observed between MC predic- 
tions and experimental data confirms a reasonable 
understanding of the background contaminations 
and reconstruction efficiency in our analysis. For 
example, after the likelihood cut, the respective 
contributions of QEL, RES, DIS and COH given 
by the cross sections and the efficiencies computed 
with the help of the MC for each process separately 
are as specified on the corresponding £ > line in 
Table m With the reduction factors for the £ > 
cut as given by the MC in the various channels, 
and normalizing the total MC to the data at this 
last stage, we find a total of 7609 events before the 
cut, distributed as shown by the figures on the cor- 
responding line. The excess of 34 MC events rela- 
tive to the data, which are necessarily mostly back- 
ground, can be taken as evidence that there is less 
than 1% excess background in the MC after the 
pre-cuts. Since the total MC background is of the 
order of 4000 events, the number found is well un- 
der the expected statistical fluctuations. Therefore, 
there is no evidence of a statistically significant dis- 
crepancy. 

In the case of 1-track events, our abilities to sup- 
press background contamination arc limited since 
all kinematic variables are expressed in terms of the 
muon momentum and emission angle 0^ with 
the help of the conservation laws for QEL events. 
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Fig. 9. The p™", a, Ou and likelihood distributions for a 
compared to real data (points with error bars). The MC distr 
data. 




mixture of QEL, RES and DIS simulated events (histograms) 
butions are normalized to the number of events observed in the 



Therefore, the proton reconstructed emission an- 
gle, Eq. can be considered as an analog of the 
likelihood function (see Fig. [T0|) . 
The explicit values for the kinematic selection cri- 
teria (£ ^ for the 2-track sample and 0.35 ^ 
dh/T^ ^ 0.5 for the 1-track sample) were found from 
the optimization of the sensitivity SG/ ^JSG + BG, 
vi^here SG and BG are the expected numbers of sig- 
nal and background events in the identified QEL 
sample. 

5.2 v^p i^i^n selection 

The investigation of antineutrino sample is a much simpler 
task since these events are mostly {^-^ 96% of cases) recon- 



structed as 1-track events (we have no hits from outgo- 
ing neutrons in the drift chambers). Therefore, we require 
identification of the positively charged muon and follow 
the procedure for the 1-track sample discussed above. The 
only difference is the absence of contamination from the 
inverse muon decay events, so we do not need to apply the 
quality cut on the transverse muon momentum. 

In Table [4] we summarize the information about the 
selection of samples with t'^n — > n~p and D^p ^^n 
candidates in the data. The last two columns of this table 
allow to make checks of compatibility between the levels 
of background in the data and in our simulations in a 
manner similar to what is explained above for the two 
track sample. 
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Fig. 10. The 6h distributions for single track i^^ (left) and Vfj, (right) samples: comparison of MC distributions (histograms) 
with the real data (points with error bars). 





57.00 GeV 




0.60 GeV^ 


w' = 


1.44 GeV^ 




0.05 GeV 








Fig. 11. A typical example of data event (run 15049 event 11514) identified as i/^n — > ^ p in this analysis. Long track is 
identified as muon, short track is assumed to be proton. 



An example of the 2-track event from real data iden- 
tified as Vfj^n — > is displayed in Fig. [TlJ 
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6 The QEL cross-section and axial mass 
measurements 

In this section we describe our analysis procedure. 
The QEL cross-section measurement using normahzation 
either to the total (DIS) (P^) CC cross-section or to 
the inverse muon decay (IMD), Vf^e^ fi^Ve, events is 
first presented in subsection 16.11 Afterwards, we describe 
the procedure used to extract the value of the axial mass 
Ma from the fit of the distribution. This is the subject 
of subsection [ 



6.1 The QEL cross-section measurement 



here N^f^ and N^ic ^^'C the number of MC events simu- 
lated and identified as QEL in the chosen detector FV. 

Similar equations can be written for any other process 
recorded in the same detector FV. If we identify A^o events 
of a process, whose flux averaged cross section in an energy 
interval containing these events is ao, we can write 

No = C^oCTo 

where <Po is the relative part of the neutrino flux belong- 
ing to the same energy interval, (we assume that iVo is 
background subtracted and efficiency corrected). 

We can now get rid of C and write the final equation 

for {<7gel)i- 



Since there was no precise knowledge of the integrated 
neutrino flux in the NOMAD experiment, we use a dif- 
ferent process with a better known cross-section, recorded 
at the same time, for the normalization of the QEL cross- 
section. A similar procedure was often applied in previous 
neutrino experiments, as for example, CERN BEBC [18] . 
Moreover, the use of another process recorded in the same 
experimental runs allows to reduce signiflcantly the sys- 
tematic uncertainty related to the detector material com- 
position. Nevertheless, this auxilliary process must meet 
two requirements: its cross-section should be measured 
with rather high accuracy and the corresponding events 
can easily be extracted from the full data sample. 

Let us divide the investigated interval of neutrino en- 
ergy into several bins and enumerate them with index 
i = 1..Ne- Then, the number of identifled QEL events 
in the i-th bin with boundaries [Ei, Ei+i] is 



Ne 



(8) 



where 



and 



= / ^{E)dE, V^, = l 

I ^^\gel{E)<P{E)dE 



Coefficient C accumulates the absolute neutrino flux 
and the number of target nucleons. The matrix element 
elf is the probability that the reconstructed neutrino en- 
ergy El, of a QEL event falls into the z-th bin, while the 
simulated energy actually belongs to the j-th bin. 

The expected background contamination is 



n'I' = C{e^{^res)+er{ad^s)) 



(9) 



where we use the deflnition of Eq. ([4]) for {(Jbg)', denotes 
the renormalized energy distribution in BG events passing 
the QEL identiflcation procedure: 



Ne 

1=1 



rec/ s 



(10) 



n: 



dat ^0^0 



A^n 



Sj ((Jdi. 



(11) 



Numerical values for {(Jres) and {cTdis) are given in Ta- 
ble [3l The efficiencies e^f , e^'^^ and ef should be esti- 
mated with the help of the MC simulation for QEL, RES 
and DIS samples separately; the factor ^oO'o/A'o comes 
from the auxilliary process used for normalization. 

Let us note that the smearing of the reconstructed neu- 
trino energy is taken into account in Eq. pip by the in- 
verse matrix of QEL efficiencies. 

Equation pT|) can also be applied to the entire energy 
interval. In this case, we can use the usual notations for 
efficiencies as in Eq. ([TU|) . From the measured {(Jqei) we 
calculate the axial mass Ma by using the Smith and Moniz 
formalism (see Fig.|4]). 

In the following subsections, we investigate the DIS 
and IMD processes which can both be used for the QEL 
cross-section normalization as just described. 

Possible sources of systematic errors in our analysis 
procedure are discussed in Section [71 



6.1.1 Selection of DIS events 

The phenomenology of neutrino DIS is well developed. Ex- 
perimental data are in rather good agreement with the- 
oretical predictions. The charged current neutrino DIS is 
an inclusive process and for its selection from the data 
sample, the following criteria are enough: 

— Fiducial volume cut. The primary vertex should be in 
the same FV as that defined for the QEL events, see 
Eq. ©; 

— Muon identification and Topology cut. At least two 
charged tracks should originate from the primary ver- 
tex; one of them should be identified as a muon (^~ 
in the case of CC and /x+ for CC); 

— Background suppression. The third criterion is used to 
avoid contributions from the QEL and RES events. We 
have checked three different possibilities for it: 
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Table 5. Selection of the DIS events in and CC samples. Total efficiency (in %), expected purity of selected events (in 
%), theoretical prediction for {adia), observed N^at and corrected A^o number of events in experimental data are given for each 
variant of DIS selection described above. 



variant of DIS selection 


1 


Vfj, sample 
2 


3 


1 


sample 
2 


3 


efficiency 


82.95 


86.84 


88.52 


75.46 


81.40 


83.20 


purity 


97.10 


98.62 


99.62 


71.48 


72.57 


73.95 


Ndat, events 


676702.0 


267517.0 


276018.0 


17744.0 


7996.0 


8500.0 


A'o, events 
relative flux #o 

(cro), 10~^** cm^ 
lO^*^ cm^ 


792162.0 
1 

16.643 
2.101 


303790.7 
0.144 
44.876 
2.127 


310617.3 
0.144 
46.069 
2.136 


16807.1 
1 

4.876 
29.012 


7128.6 
0.106 
20.124 
29.924 


7553.4 
0.106 
21.999 
30.872 



1. The total visible energy in the event should be 
El, ^ 300 GeV and the reconstructed hadronic 
mass W ^ lA GeV; in this case the computation 
of {adts) has been done for GRV98-LO PDF model 
according to the prescriptions in [57] . 

2. We keep the requirement for the reconstructed had- 
ronic mass (W ^1.4 GeV) but reduce the neutrino 
energy region to 40 ^ i?^ ^ 200 GeV; theoreti- 
cal calculation of (adis) is also done with the help 
of [57]. 

3. Using the same neutrino energy interval as in 2. 
(40 ^ ^ 200 GeV), we remove the cut on the 
reconstructed hadronic mass W. In this case, we 
take the total CC neutrino-nucleon cross-section to 
be:0 

al°\E^)lE^ = (0.677± 0.014) x 10"^* cmVGeV 
ctI°\E^)/E^ = (0.334 ±0.008) x 10"^* cmVGeV 

(PDG average [49]). The calculated {a tot) should 
be corrected due to the fact that NOMAD target 
is slightly non-isoscalar. 

The numerical results of the DIS events selection can 
be found in Table For the QEL cross-section normaliza- 
tion we use results obtained with the last method (PDG 
based) as having the most solid ground. Thus, the final 
normalization is performed to the total (j?^) CC cross- 
section. We also checked that this normalization is con- 
sistent with two previous calculations based on approach 
from [57] within 1.6% (5.9%) for j/^ (i?^) CC sample. 

6.1.2 Selection of inverse muon decay events 

Inverse muon decay fc'^e^ li^Ve is a purely leptonic 
process, which is well known both on theoretical and ex- 
perimental grounds. Its cross-section in the Born approx- 
imation is: 

/ < V 

<J«nd{E,) = CTasE, 1 - — ^ (12) 



^ The CHORUS measurement for the CH2 target [81] is con- 
sistent with this value. 



The numerical value of the constant CTos calculated in the 
framework of the Standard Model was found to be in good 
agreement with experimental measurements f82| : 

^ fUl^ ^ 1.723 X 10-41 cm2 GcV-^ (13) 

TT 

The number of IMD events Nq is proportional to its 
flux averaged cross-section from Eq. (|4]) : 

(aw) = 1.017 X 10-"*" cm2 (14) 

and expected to be at least 650 times smaller than the 
number of DIS events. 

To select the IMD events we require: 

— the primary vertex should be in the same fiducial vol- 
ume as that used for identified QEL events, sec Eq. ©; 

— there is only one negatively charged track originating 
from the primary vertex; it should be identified as a 
muon; 

— there are no veto chamber hits in the vicinity of the 
intersection point of the extrapolated muon track and 
the first drift chamber (quality cut, the same as for 
1-track events from the QEL sample); 

— the muon energy is above the threshold: 

En > — = 10.93 GeV (15) 

2me 

— the transverse momentum of the muon produced in 
IMD event is very limited by kinematics: 

In this sample the contamination from the reaction 
DeS ^J■^^^fJ, is estimated to be at the level of ^ 10^^, e.g. 
well below 1 event, since the ratio of the fluxes i^e/^'p is 
0.0027 [SO] while the ratio of the cross-sections is a{Dee — > 
H~Pf^)/a{Vf_,e H~Ve) ~ 1/3. 

We determine the number of signal events Ni^d from 
the fit of the p\ distribution to experimental data with 
the function F{p\): 

F{pI) = N,„,dF^^d{p\) + [Ndat - N,„,d] Fbgipl) (16) 

where Fimd and Fbg are the normalized MC expectations 
for signal and background p"^ distributions; Ndat denotes 
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• NOMAD data 

1 parameters fit 

^ IMD signal 




pi, (GeV/c)^ 



Fig. 12. Inverse Muon Decay: NOMAD experimental data, 
the non-equidistant p\ distribution. 



the number of events in real data which passed aU selection 
criteria. 

The QEL events are now playing the role of the most 
important background for the IMD selection. However, 
the contaminations from the RES and DIS events cannot 
be neglected since they distort the shape of the distri- 
bution. As usual, the relative contribution of each process 
to the expected background is proportional to the cor- 
responding efficiency and flux averaged cross-section (see 
Table O. 

The expression p6p contains only one free parameter 



N,, 



which is the number of observed IMD events. Fi- 



nally, for p-]_ < 0.03 (GeV/c)^ interval we find iVimd = 
436.0 ± 28.5 with the quality of the fit x^/NDF = 0.89 
(see Fig. [T2|) . Taking into account that the selection ef- 
ficiency for the IMD events is 87.8% we report the total 
number of IMD events Nq, which can be used for the QEL 
normalization: 



iVn = 496.6 ± 32.5 



(17) 



6.2 Axial mass measurement from the distribution 

To extract the axial mass from the distribution the 
experimental data are fitted to the theoretical predictions 
using a standard method. We bin the events in two 
variables and E^, (in the case of a single interval our 
procedure can be considered as the usual 1-dimensional 
fit)0. 

Let us enumerate bins with index i ~ 1..Nb', bin i ~ 
Nb + 1 contains events which fall outside of the investi- 
gated {Efj,Q^) region. It is convenient to define bound- 
aries in such a way that each bin with i ~ I..Nb contains 
approximately the same number of experimental events 
passing all identification criteria. 

A minimization functional is: 



dat 



(18) 



where Nf""*^ is the number of events in the i-th bin of 
the non-weighted experimental distribution, while A^/'' is 
a superposition of the normalized MC background A^^^^ 
and the expected QEL signal: 



A^, 



th 



{Ma) = A^^ 



N^^ being defined 



(19) 



is the probability that 



This equation is similar to Eq. 
in the same way as in Eq. ([5]); e^j' 
a QEL event simulated in the j-th bin is reconstructed in 
the i-th bin. The QEL scattering dynamics is described 
by the following term: 



1 



da 



{E,Q^,MA)<P{E)dEdQ^ (20) 



qel )i I 



-Ne 



Nb 
J = l 



(21) 



here Hi denotes the {E^,Q^) interval, which corresponds 
to the i-th bin; da/dQ^ is the differential QEL cross- 
section on bound target nucleon (see formulae in [32]) 
The coefficient C can be defined in either of two ways: 

1 . the A*'* distribution is normalized to the total number 
of events in the experimental data: 



The relative error for (Tq/Nq in the IMD case is about 
7% (due to the small statistics of the IMD sample). Nev- 
ertheless the normalization factor itself, = 2.048 x 
10~*^ cm^, is in agreement (within ^ 4%) with the eval- 
uation based on the DIS sample (see Table [5]). 

The use of the IMD process for the normalization is 
an interesting independant cross-check of the absence of 
problems in our procedure. In particular, it allows to ver- 
ify that there are no effects arising from possible trigger 
inefficiencies in the selection of neutrino events consisting 
of a single muon going through the trigger planes. 



Nb 



Nb 
i=l 



(22) 



In this case, the proposed method should be sensitive 
only to the shape of the distribution but not to the 
absolute number of identified events (contrary to the 
Ma measurement from the total QEL cross-section). 



* In practice it is convenient to use dimensionless variables 
(a, h) instead of {E^, Q^). Then, = E^'" + a{ElJ"'=' - E^'") 
and Q2 = Ql^i„{E,) + b[Ql,^^{E,) - Ql,,„{E,)]. So, a,b € 
[0,1]. 
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Fig. 13. The Q distributions in identified QEL events. 

2. C is defined in the same way as for the total QEL 
cross-section measurement, i.e. we use another process 
(DIS) for normalization: 



C 



No 



(23) 



If we sum over the variable for the investigated 
(-Ey, Q^) interval, finding the Ma parameter from Eq. (fTSi 
becomes nothing else than the numerical resolution of 
Eq. ([8|). Therefore, this variant of the fit can be con- 
sidered as a simultaneous fit of the total and differen- 
tial cross-sections; henceforth, we shall refer to it as 
a (g) da/dQ'^ fit. 

Fig. [13] presents a comparison of the reconstructed 
distribution with our MC prediction. The expected back- 
ground contamination is also shown. 

We can now apply the proposed methods to experi- 
mental data and measure the QEL cross-section and axial 
mass Ma - The numerical results are reported in Section[8l 
while the discussion of the corresponding uncertainties is 
presented in the next section. 



7 Systematic uncertainties 

We have studied several sources of systematic uncertain- 
ties, which are important for the measurement of the to- 
tal QEL cross-section and axial mass parameter. They are 
hsted below: 

1. identification of QEL events; we vary the selection cri- 
teria within reasonable limits (£ > ± 0.4 for 2-track 
sample and Bii/tt > 0.35 ± 0.03 for 1-track sample). 
The final result is found to be practically insensitive to 
the exact positions of the muon azimuth </j^ cut and ad- 
ditional requirements for the p™^ ^ a and Oh variables: 



e.g. in the i/^ analysis a more strict cut O.Itt < i/j^ < 
O.Ott leads to 0.8% variation in the measured cross sec- 
tion while a change in the prc-cuts to P™** < 0.9 GeV, 
a/iT > 0.75 and O.IStt < 6^ leads to an uncertainty of 
0.4%. 

uncertainty in the total (mainly DIS) charged current 
muon neutrino cross-section, which enters both in the 
normalization factor cfo/Nq and in the subtraction of 
the corresponding DIS background (the experimental 
error on (adis) is 2.1% for i/^ CC and 2.4% for CC); 
uncertainty in the RES cross section, which determines 
the contamination admixture of the single resonant 
pion events in the identified QEL sample (we assume 
10% error on {(Tres) both for neutrino and antineutrino 
cases, see e.g. [56]): 

FSI interactions (we vary tq and a^^^ DPMJET pa- 
rameters for fixed A/^"= = 1.03 GeV); 
uncertainty in the neutrino flux shape (the relative er- 
rors for each bin were taken from [3D]); 
neutral current admixture (we assume 5% error for 
the corresponding cross section, which can be foimd in 
Table [3]); 

charge misidcntification of the primary Icpton (recon- 
structed i^^ CC event is classified as Df^ CC and vice- 
versa) ; 

contamination from coherent pion production (see sub- 
section 14.41) . 



In Table [6] we present our numerical estimations for 
systematic uncertainties (in the case of scattering, sys- 
tematic errors were calculated for the mixture of 1-track 
and 2-track subsamples). One can sec that the most im- 
portant contributions come from the QEL identification 
procedure and from the uncertainty on the non-QEL pro- 
cesses contribution to the selected sample of signal events. 

The nuclear reinteractions (FSI effect) significantly af- 
fect the neutrino sample only (see Table [5]), while in the 
antineutrino case the infiucnce of the nuclear reinterac- 
tions is expected to be negligible. For scattering, the 
cross-sections can be calculated separately for both the 1- 
track and 2-track subsamples of identified QEL events or 
for their mixture. We can then compare the results and 
choose whichever one has the minimal total error. In our 
case it was obtained for the combined 1-track and 2-track 
sample, which was found to be almost insensitive to the 
variation of DPMJET parameters (see Section [8] for ex- 
planations). 

The uncertainty on the shape of the (anti)neutrino 
spectrum is important for the measurement of <7qei as a 
function of neutrino energy E^, . But it does not affect both 
the flux averaged cross section {<Jqei) and the Ma extrac- 
tion from the distribution. 

The uncertainty due to the primary lepton misidenti- 
fication and neutral currents comes into play through the 
subtraction of the corresponding background from the se- 
lected DIS sample, that is, from the normalization factor. 
The admixture of those events into the identified QEL 
events is negligible. 
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Table 6. The relative systematic uncertainties (in %) of the QEL cross section {(Jqei) and axial mass Ma, measured in 
i^fin fi^p and Ufj,p — » fi'^n reactions. 





Source 




Ma from (cr,ei>^„ 


Ma from dayjiiQ^ 




Ma from {(Jqd)ii^ 


1 


QEL identification procedure: 














likelihood or 6^ cut 


3.5 


3.2 


2.4 


4.3 


4.2 




<Pi_i cut 


0.8 


0.7 


0.3 








P™", Q and 9h precuts 


0.4 


0.4 


0.4 






2 


5{(Jdis) 


2.9 


2.6 


0.2 


4.2 


4.2 


3 


5(o"res) 


4.0 


3.6 


0.6 


7.6 


7.4 


4 


nuclear reinteractions 


1.8 


1.6 


6.5 






5 


shape of neutrino spectrum 


0.2 


0.2 


0.1 


0.9 


0.9 


6 


NC contribution 


< 0.1 


< 0.1 




1.1 


1.1 


7 


muon misidentification 


< 0.1 


< 0.1 




1.0 


1.0 


8 


coherent pion production 


< 0.1 


< 0.1 


< 0.1 


1.1 


1.1 




total 


6.5 


5.9 


7.0 


9.9 


9.5 



1.4 

M 

o 1 

D 0.8 
0.6 
0.4 
0.2 




3 NOMAD 1 -track sample 
• NOMAD 2-track sample 
■ NOMAD combined sample 
- □ M, error 0.06 GeV 



\^-\-n-> \i +p 




10 



10 



(GeV) 



Fig. 14. Comparison of our {aqei)v^ measurements as a function of the neutrino energy in the 1-track and 2-track subsamples 
(for the best parameter ro = 1.0) with the final {ag^i}^^ values measured using the full event sample, see Table [7] 



8 Results 

8.1 i^f^n fjL^p sample 

The results of our analysis for the sample are sum- 
marized in Table [9l We measure the flux averaged QEL 
cross-section in the neutrino energy interval 3 — 100 GeV 
(see Eq. ([TIT) ) for the 1-track and 2-track samples as well 
as for their mixture (which is called Combined in Table [9]). 
For each {(Jqei) we calculate the corresponding axial mass 
value, Ma. Results on Ma extraction both from the stan- 
dard fit and from the combined a®da/dQ^ fit arc also 
given. These measurements are repeated for several QEL 
MC with different values of input parameters (the axial 
mass Ma was varied between 0.83 and 1.23 GeV in steps 
of 0.1 GeV; the formation time tq was allowed to take a 
value of 0.6, 1.0 and 2.0; the correction factor a^g^ was 
varied within the interval [0.54,0.69]). On top of this the 
NUANCE QEL MG with its own treatment of FSI effects 
is used for cross-checks. 

We then observe that Ma recalculated from the mea- 
sured {(Tqei) depends on tq if one refers to the 1-track or 
the 2-track samples. Specifically, the measured Ma value 



increases with increasing tq when extracted from the 1- 
track sample while it decreases when extracted from the 
2-track sample. This can be understood if we take into 
account the fact that the tq parameter controls the proba- 
bility for an outgoing nucleon to be involved in an intranu- 
clear cascade. Increasing tq then increases the fraction of 
QEL events with reconstructed proton and thus populates 
the 2-track sample to the detriment of the 1-track sample. 
This is the reason for the systematic overestimation of Ma 
extracted from the 1-track sample alone and its underes- 
timation when extracted from the 2-track sample alone. 
However the value of Ma extracted from the combination 
of the 1-track and 2-track samples is almost insensitive to 
variations of the tq parameter. 

We also find that using the QEL Monte Garlo with 
To = 1 and Qf^Q^ = 0.6 provides the most accurate predic- 
tion for the ratio between the 1-track and 2-track samples 
(and hence the most adequate description of the FSI): in 
this case the flux averaged QEL cross-section stays ap- 
proximately the same whether measured from the 1-track 
sample or from the 2-track sample (see Table O . This al- 
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Fig. 15. Comparison of NOMAD (o-q^i)^^ measurements witii previous experimental data on Vf^ scattering off lieavy nuclei 
(ANL 69 (Spark-chamber) [1], NuTeV 04 (FermiLab) [23], CERN HLBC 69 (CERN, Heavy Liquid Bubble Chamber) [14], 
CERN GGM 77 (CERN, Gargamelle BC) [H], CERN GGM 79 [13, IHEP 85 (IHEP, spark-chamber) [20], IHEP SCAT 90 
(IHEP, BC) [22]). The solid line corresponds to the Ma value obtained in the NOMAD experiment, the error band takes 
into account both statistical and systematic uncertainties of the present analysis. Nuclear effects are included into calculations 
according to the relativistic Fermi gas model by Smith and Moniz |71j for Carbon with binding energy Et = 25.6 MeV and 
Fermi momentum Pp — 221 MeV/c. 
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Fig. 16. Comparison of NOMAD {oqei)v^ measurements with previous Vfi,D experimental data (ANL 73 (Argonne 12- foot 
BC) [2], ANL 77 [3], BNL 81 (Brookhaven 7-foot BC) [6], FNAL 83 (FermiLab 15-foot BC) [9], BEBC 90 (CERN, Big European 
Bubble Chamber) [18]; corrections for nuclear effects have been made by the authors of the experiments). The solid line and 
error band corresponds to the Ma value obtained in the NOMAD experiment. 



lows US to exclude the MC sets with tq = 0.6 and 2.0 from 
further considerations. 

Fig. [T4l shows a comparison of om' {aqei)^^^ measure- 
ments as a function of the neutrino energy in the 1-track 
and 2-track subsamplcs (for the best parameter tq = 1.0) 
with the final {cTqei)^^ values measured using the full event 
sample. 

Similarly we have observed that when using the full 
sample (1-track and 2-track) the measured Ma is not very 
sensitive to modifications of the a^^^ parameter. And us- 



ing the NUANCE simulation code as a cross check gives 
a very consistent picture: the Ma value extracted from 
the 1-track sample is also different from the one extracted 
from the 2-track sample, while the value obtained with 
the combined sample nicely agrees with our measurement 
with the best FSI parameters. Thus, our results for the 
neutrino case are: 

{aqei)v,, = (0.92±0.02(siai)±0.06(syst)) x 10"^^ cm^ 
Ma = 1.05 ± 0.02(stat) ± 0.06(syst) GeV (24) 
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Fig. 17. Comparison of NOMAD {aqsi)v^ measurements, Table |8l witii previous experimental data on i/^ scattering off heavy 
nuclei (CERN GGM 77 CERN GGM 79 [1^, IHEP 85 [20], IHEP SCAT 90 ^ and NuTeV 04 [23]). The solid line and 
error band corresponds to the Ma value obtained in the NOMAD experiment, the error band takes into account both statistical 
and systematic uncertainties of the present analysis. Nuclear effects are included into calculations according to the standard 
relativistic Fermi gas model. 



This result ([^^ is indeed in agreement with both the stan- 
dard fit of the distribution: 

Ma = 1.07 ± 0.06{stat) ± 0.07{syst) GeV (25) 

and the fit of the combined a^da/dQ"^ distribution of the 
NOMAD data: 

Ma = 1.06 ± 0.02(siai) ± 0.06(syst) GeV (26) 

(see Table [HI these results are obtained with a QEL MC 
using Ma = 1.03 GeV). 

We use the 2-track sample only to extract Ma from 
the fit of the distribution since in this case the purity 
of QEL identification is rather high 74%, see Table [?]). 
The results depend on the input MC parameters (axial 
mass and formation time) but still are in nice agreement 
with the results of the extraction of Ma from the measured 
QEL cross-section based also on a 2-track sample analysis. 
This can be considered as an additional confidence for our 
measurements using the full QEL sample. 

The measured cross-section of the f^n — > ^,~p reac- 
tion as a function of the neutrino energy is presented in 
Table [7] and is shown in Figs. [TSl and [TBI These results are 
compared to the previous measurements performed with 
deuterium and heavy nuclei targets (see discussion in Sec- 
tion [2]). 



8.2 v^p fi^n sample 

In the case the event topology is just a single /i+ , thus 
the uncertainties in the treatment of FSI effect almost do 
not influence the event selection. Since our measurement 
of the cross-section of the Df^p jX^n reaction is based 
on a 1-track sample only, we do not show the dependence 



of the results on the variation of the tq and a^^^ pa- 
rameters. Instead we display a dependence on the input 
Ma in Table [TO] The results for the measured Ma are 
found to be quite stable. In Fig. [17] we show the mea- 
sured Vfj,p fi'^n cross section as a function of the an- 
tineutrino energy superimposed with the theoretical curve 
drawn with Ma = 1.06±0.12 GeV and with nuclear effects 
according to the standard relativistic Fermi gas model. 
Table [5] summarizes our results for the D^p —f fi'^n cross- 
section measurement in the different antineutrino energy 
intervals. The cross-sections are measured on a Carbon 
target and also recalculated for a free nucleon. The sta- 
tistical and systematic errors are both provided. The ob- 
served number of events in the data, the predicted number 
of background events, the background subtracted and ef- 
ficiency corrected number of events are also shown. 
Our final results for the antineutrino case are; 

{(Jqei)i?^ = (0.81 ± 0.05{stat) ± 0.08{syst)) x 10"^^ cm^ 
Ma = 1.06 ± 0.07(stat) ± 0.10{syst) GeV (27) 

9 Conclusions 

The cross-section measurement of the z/^n — > n~p and 
D^p fi^n reactions on nuclear target was performed and 
reported in this article. The samples used in the analysis 
consist of 14021 neutrino and 2237 antineutrino events, 
which were identified as quasi-elastic neutrino scattering 
among the experimental data collected by the NOMAD 
collaboration. 

We have discussed in details the analysis procedure 
and the most significant sources of systematic error. Spe- 
cial attention was paid to the infiuence of the FSI effects 
on the measured physical quantities. The DPMJET code 
was used to simulate these FSI effects. We also proposed 
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Table 7. Cross-section of quasi-elastic neutrino scattering (in units of 10^^* cm'^, statistical and systematic errors), acarton 
is measured for NOMAD nuclear target and normalized per 1 neutron; a Free ~ cTNuci/g is the cross-section for the free target 
neutron (the factor g is calculated according to the Smith-Moniz model, see [32]). The number of selected events in raw data 
Ndat, the estimated background contamination Ntg and the number of events Near corrected for background and efficiency are 
also reported. The difference in the total number of data events with respect to Table |4] (13683 vs 14021) is due to the additional 
cut on the neutrino energy 3 < Ev{GeV) < 100. 









Ndat 


Nt, 


cor 


(a ± Sa^tc 


It ± Sasy 


st)Nucl 


{a ± Sastat 


± Sasy 


st^Free 


3 - 


6 


4.7 


396 


211.6 


660.6 


0.994 


0.125 


0.078 


1.057 


0.133 


0.083 


6 - 


9 


7.7 


1115 


580.2 


1663.4 


0.942 


0.072 


0.078 


1.001 


0.077 


0.083 


9 - 


12 


10.5 


1683 


835.3 


2591.1 


1.014 


0.061 


0.075 


1.077 


0.065 


0.080 


12 - 


15 


13.5 


1647 


834.9 


2310.7 


0.859 


0.057 


0.070 


0.913 


0.060 


0.075 


15 - 


21 


17.8 


2815 


1451.6 


3766.8 


0.843 


0.040 


0.067 


0.896 


0.043 


0.071 


21 - 


27 


23.8 


2040 


956.2 


3084.7 


0.991 


0.052 


0.070 


1.053 


0.055 


0.075 


27 - 


33 


29.8 


1279 


610.5 


1816.8 


0.888 


0.064 


0.073 


0.943 


0.068 


0.077 


33 - 


39 


35.8 


852 


400.9 


1246.2 


0.970 


0.086 


0.078 


1.030 


0.091 


0.083 


39 - 


54 


45.3 


1008 


496.1 


1397.8 


0.896 


0.070 


0.074 


0.951 


0.074 


0.078 


54 - 


100 


71.7 


848 


416.9 


1191.5 


0.911 


0.073 


0.077 


0.967 


0.078 


0.082 


3 - 


100 


23.4 


13683 


6794.2 


19718.2 


0.919 


0.017 


0.060 


0.976 


0.018 


0.063 



Table 8. The same as Table [7] but for antineutrino. 









Ndat 


Nby 


cor 


(cr ± Sastc 


It zb Sasyst) Nucl 


(a ± Sastc 


It ± Sasy 


st'}Free 


3 - 


6 


4.5 


291 


176.4 


181.9 


0.585 


0.097 


0.079 


0.639 


0.106 


0.086 


6 - 


9 


7.4 


292 


182.6 


159.5 


0.710 


0.125 


0.093 


0.767 


0.135 


0.100 


9 - 


15 


11.8 


464 


286.2 


254.1 


0.851 


0.114 


0.096 


0.913 


0.122 


0.103 


15 - 


24 


19.0 


425 


274.0 


210.8 


0.822 


0.125 


0.095 


0.878 


0.134 


0.102 


24 - 


42 


31.5 


411 


252.5 


223.0 


1.009 


0.141 


0.098 


1.075 


0.150 


0.104 


42 - 


100 


60.1 


288 


181.9 


145.9 


0.993 


0.173 


0.099 


1.056 


0.184 


0.105 


3 - 


100 


19.7 


2171 


1353.6 


1182.5 


0.811 


0.053 


0.081 


0.866 


0.056 


0.086 



a method for tuning the intranuclear cascade parameters 
(mainly the formation time tq), which was then used to 
reduce the corresponding systematic uncertainty. 

For the z^^ case stable results have been obtained with 
the combined 1-track and 2-track samples since they are 
almost insensitive to the FSI effects. 

The results for the flux averaged QEL cross-sections 
in the (anti)neutrino energy interval 3 — 100 GeV are 
Wqei)^^ ^ (0.92±0.G2(stat)±G.06(syst)) x IQ-^s cm^ and 
{agei}i?^ = (0.81±0.05(stat)±G.08(syst)) x IQ-^s cm^ for 
neutrino and antineutrino, respectively. 

The axial mass Ma was calculated from the measured 
cross-sections: we find Ma = 1.05 ±0.06 GeV from the i/^ 
sample and Ma = 1-06 ± 0.12 GeV from the sample. 
The Ma parameter was also extracted from the fit of the 

distribution in the high purity sample of i'^ quasi- 
elastic 2-track events (with a reconstructed proton track). 
It was found to be consistent with the values calculated 
from the cross-sections. 

Our results are in agreement with the existing world 
average value [331136] and do not support the results found 
in recent measurements from the NuTeV [23], K2K [24l 
[25] and MiniBooNE [26] collaborations, which reported 
somewhat larger values, however still compatible with our 
results within their large errors. A summary of existing 



experimental data on the axial mass measurements in neu- 
trino and antineutrino experiments is shown in Fig. 1181 

It should also be noted that the preliminary results 
reported earlier by the NOMAD collaboration for the 2- 
track sample only [5^[S3] suffered from a large systematic 
bias related to an improper treatment of the FSI effects in 
the simulation program. They should be now superseeded 
by the new measurements reported here. 
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Fig. 18. A summary of existing experimental data: tlie axial mass AIa as measured in neutrino (left) and antineutrino 
(right) experiments. Points show results obtained both from deuterium filled BC (squares) and from heavy liquid BC and other 
experiments (circles). Dashed line corresponds to the so-called world average value A'Ia ~ 1.026 ± 0.021 GeV (see review [33] 1. 
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Table 9. Parameters of the QEL MC simulation (axial mass M^'^ and parameters of FSI modeling) are listed in the first three 
columns. The intermediate columns contain results of the QEL cross-section measurement (in units of 10~^*cm^, without 
errors) for the different topology of identified events (with or without reconstructed proton track). The axial mass value obtained 
from the fit of distribution and a ® da/dQ^ fit {Qfi^ — 0.2 — 4GeV^) are given in the last columns of the table; the given 
statistical errors on the axial mass 5Ma are from MINUIT output. 



To 


MC parameters 

^ F n ,rmc 

amod Ma 


Single 

O'qel 


track 
Ma 


Two tracks 

Cqel Ma 


Co 

Cqel 


mbined 
Ma 


Fit of distribution 
Ma ± SMa 


Fit of a® da/dQ^ 
Ma ± SA4a 


0.60 
1.00 
2.00 
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0.60 


0.83 
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0.83 
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0.885 
0.918 
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1.015 
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1.014 
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0.977 
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0.057 
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19.3 
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0.60 


0.93 
0.93 
0.93 
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0.931 
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0.844 
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0.968 
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0.896 
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1.026 


1.135 
1.075 
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0.018 
0.018 


27.1 
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1.092 
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0.077 
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0.017 


26.4 
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0.60 
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1.081 


0.926 
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0.940 


1.073 


1.092 


0.056 
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18.2 
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0.994 


1.127 
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1.088 
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0.60 
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1.035 
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Table 10. The results of QEL cross-section measurement. The parameters of the DPMJET model are tq = 1.0, cJmod = 0.6. 
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